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We study the phase diagram of the asymmetric Hubbard model (AHM), which is characterized 
by different values of the hopping for the two spin projections of a fermion or equivalently, two 
different orbitals. This model is expected to provide a good description of a mass-imbalanced cold 
fermionic mixture in a 3D optical lattice. We use the dynamical mean field theory to study various 
physical properties of this system. In particular, we show how orbital-selective physics, observed 
in multi-orbital strongly correlated electron systems, can be realized in such a simple model. We 
find that the density distribution is a good probe of this orbital selective crossover from a Fermi 
liquid to a non- Fermi liquid state. Below an ordering temperature To, which is a function of both 
the interaction and hopping asymmetry, the system exhibits staggered long range orbital order. 
Apart from the special case of the symmetric limit, i.e., Hubbard model, where there is no hopping 
asymmetry, this orbital order is accompanied by a true charge density wave order for all values of the 
hopping asymmetry. We calculate the order parameters and various physical quantities including 
the thermodynamics in both the ordered and disordered phases. We find that the formation of the 
charge density wave is signaled by an abrupt increase in the sublattice double occupancies. Finally, 
we propose a new method, entropic chromatography, for cooling fermionic atoms in optical lattices, 
by exploiting the properties of the AHM. To establish this cooling strategy on a firmer basis, we also 
discuss the variations in temperature induced by the adiabatic tuning of interactions and hopping 
parameters. 
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I. INTRODUCTION 

Strongly correlated electronic systems are known to ex- 
hibit a wide variety of interesting and novel phases rang- 
ing from Fermi liquids, heavy fermions, Mott insulators, 
superconductors to other states with complex long range 
orders. The field of strongly correlated fermions in con- 
densed matter shares many common research goals with 
that of cold atoms loaded in optical lattices. Yet these 
two fields emerge from different scientific communities 
and often look at similar problems with different tools 
and from different physical perspectives. In this context, 
the present theoretical study focuses on a model hamil- 
tonian that generalizes central models of strongly corre- 
lated electron systems, which, nevertheless, might only 
be experimentally realized with cold atoms in an optical 
lattice. The present work is a contribution towards build- 
ing bridges between these two scientific communities. 

The paradigmatic model of correlated fermions is the 
Hubbard model^ (HM), which describes itinerant elec- 
trons on a lattice with hopping t and subjected to an 
onsite repulsion U stemming from screened Coulomb in- 
teractions between two electrons. This model, though 
the focal point of interest in condensed matter physics 
for decades, however remains unsolved in two and three 
spatial dimensions. In the limit of infinite dimensions, 
or infinite lattice connectivity, the model remains non- 
trivial but is exactly solvable using the dynamical mean 



field theory (DMFT)^. DMFT shows that in the absence 
of any frustration, the ground state of the system at half- 
filling is an antiferromagnetic insulator for all values of U 
with an associated Neel temperature. In cases where due 
to inherent frustration the system is in a paramagnetic 
state, one finds a correlated metal described by a Fermi 
liquid with heavy quasiparticles at low and moderate U . 
Upon increasing the interaction strength, always at half- 
filling, the heavy metal breaks down and a Mott insulat- 
ing state is realized with a coexistence region^. Such a 
Mott transition was indeed observed in the classic three 
dimensional Mott-Hubbard system V203^. 

Another important model Hamiltonian of strongly cor- 
related systems is the Falicov-Kimball modeP (FKM), 
which is a simplified variant of the HM where one of 
the spin species has zero hopping, hence, its motion is 
frozen. This model explicitly breaks the spin SU(2) sym- 
metry and DMFT studies^ showed that it has an anti- 
ferromagnetically ordered insulating phase at half-filling, 
and regions of phase separation at low values of the onsite 
Coulomb interaction for finite doping. Similar to the HM, 
a paramagnetic Mott metal -insulator transition (without 
any coexistence) as a function of the interaction strength 
was also found in the FKM at half filling. A fundamen- 
tal difference between the two models is that while the 
metallic state in the HM was described by a Fermi liq- 
uid, the metallic state in the FKM is a non-Fermi liquid 
with no quasiparticles. This exotic state stems from the 
broken translational symmetry due to the random distri- 
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buttons of the immobile electrons on the lattice. 

The natural connection between these two models is 
the Asymmetric Hubbard model (AHM), where each spin 
species has a different hopping amplitude. This model, 
originally considered purely academic in solid state sys- 
tems, is of particular relevance to cold atoms on optical 
lattices. In fact, asymmetric Hubbard models arc realiz- 
able in optical lattices, using either two hyperfine states 
of a spin- 1/2 atom or two different spinless fermionic 
atoms (e.g., ^Li and ^°K), as described in Refs. 6-12. 
Each of the hopping amplitudes ta can be independently 
controlled by changing the intensity of the lasers and the 
Coulomb interaction strength U is related to the scatter- 
ing length and can be well controlled by Fano-Feshbach 
resonance^'^'^*. Following Refs. 15 and 16, the hopping 
asymmetry r in a mixture of ^Li/^°K can be tuned from 
values close to the FKM (r=0) to that of the HM (r=l). 

Another important reason to study the AHM, is the 
fact that it can be thought of as a minimal model of multi- 
orbital systems. Earlier multi-orbital models in (bosonic) 
cold-atom systems have been considered in Refs. 17-19. 
In general, in correlated materials, often when there is 
more than one band close to the Fermi energy, new phys- 
ical phenomena arise. The existence of many interacting 
bands close to the Fermi energy, may generate exotic 
phases, produced by the difference in the population of 
each band, the inter and intra band Coulomb interaction, 
as well as the Hund's coupling. Examples of these mate- 
rials, which are subjects of current research include high- 
Tc iron-based superconductors, where 5 bands cross the 
Fermi energy. From this perspective, the AHM, which 
can be viewed as a model for spinless fcrmions with two 
different orbitals, can be considered as an attempt to ap- 
proach the physics of correlated multi-orbital systems, in 
cold-atoms. For this reason, in this paper we adopt the 
notation of two different orbitals h and /, instead of two 
spin-projections f and J,, as we will explain later. 

Most of the recent theoretical work on the AHM in 
the cold-atoms context, has concentrated on establish- 
ing the phase diagram in the one dimensional case, for 
both, attractive and repulsive interactions. A rich va- 
riety of groundstates have been found, including Mott 
insulator, charge density wave and superconductivity, in 
the case of mixtures of equal density^, and also FFLO 
states in the case of population imbalance^°. Beyond the 
one dimensional case, the AHM with attractive interac- 
tions on a cubic optical lattice, has been studied using 
DMFT^^. A related, but more general model, was also 
analyzed in Rcf. 20 using the slave-boson technique. The 
phase diagram of the repulsive AHM without any long- 
range order was studied using DMFT in Refs. 16 and 21. 
Both studies found a Mott transition with a region of two 
coexistent solutions, at low enough temperatures for all 
non-zero values of the hopping asymmetry. In particu- 
lar, in Ref.21, we showed that an interesting temperature 
driven orbital-selective crossover takes place in the AHM. 
At temperatures below a coherence temperature Tcoh, the 
system is a heavy Fermi-liquid, qualitatively analogous to 



the Fermi liquid state of the symmetric Hubbard model. 
Above Tcoh, instead of the standard incoherent metal 
seen in the HM case, an orbital-selective crossover oc- 
curs, wherein one fermionic species effectively localizes, 
and the metallic state maps onto the non-Fermi liquid 
of the Falicov-Kimball model. This is a validation of 
the fact that the AHM can be considered as a minimal 
model for the study of orbital physics. Various observ- 
ables such as the double occupation, the specific heat and 
entropy were also computed. As stated earlier, a phase 
diagram without any long-range order is appropriate pro- 
vided an inherent frustration prevents the formation of 
any kind of orbital ordering. However, for the unfrus- 
trated case, which corresponds to the bare model hamil- 
tonian without additional assumptions, one needs to ex- 
plore ordered solutions to the mean field equations. Here, 
for simplicity, we restrict ourselves to the case of half- 
filling, where only Neel-type (ie checkerboard) order is re- 
alized. However, away from half-filling, incommensurate 
order may also occur^'^'^. To explore the long-range or- 
dered solutions and to understand the way how these are 
connected to the previously studied finite temperature 
crossover regimes, including the orbital selective state^^, 
is the main goal of the present work. In this paper, we 
study the orbital-ordered phases as a function of hopping 
asymmetry, interaction strength and temperature, which 
then permits us to obtain a complete phase diagram of 
the half-filled AHM. Since fermionic cold atom systems 
have not yet attained the low temperatures which are 
requisite for observing any kind of order, we also propose 
a novel cooling scheme which is based on the hopping 
asymmetry of the two orbitals. 



The present problem is at the interface between two 
different fields, thus, we summarize all the technical de- 
tails in a self-contained section II A that can be skipped if 
desired. For further information on the DMFT-technique 
we suggest Ref. 2. Qualitative features of the long-range 
order phase are described in II B. In section HI A, we 
briefly review the results of the disordered phase of the 
AHM found in Ref. 21, and discuss their consequences 
for cold atoms experiments. In section IIIB, we discuss 
the ordered phase of the AHM and present results for 
the critical temperature To , the different order param- 
eters and the thermodynamical quantities of the model. 
In section HI C, we propose a strategy for cooling down a 
mixture of cold fermionic atoms, that exploits the mass 
imbalance and may have practical applications. In sec- 
tion HID, we analyze the effects on the temperature of 
the condensate of tuning the values of hopping and in- 
teraction of the atomic particles. Finally in section IV, 
we summarize our conclusions. The expressions used to 
compute the physical observables can be found in the 
appendix. 
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II. MODEL AND METHOD 
A. Mean-field equations 

The asymmetric Hubbard model Hamiltonian reads, 
H = - tAaCja + uY.ina^l/2){n,u~l/2) 

<i.j>](T—l,h i 

(1) 

where cj^. , Cia- are the fermion creation and annihilation 
operators at site i and a = l,h labels the two orbitals, 
light and heavy. t„ and U denote the hopping and the re- 
pulsive interaction strengths and the asymmetry param- 
eter is defined as r=th/ti. The limits r=0 and r=l cor- 
respond to the FKM and HM, respectively. We restrict 
our study to the particle-hole symmetric half-filled case, 
with number of particles per site (n;) = (n/i) = 1/2. To 
describe a realization of the AHM in cold atom systems 
in an optical lattice, one would need to add a confining 
potential to the hamiltonian (1). For simplicity, in this 
initial study, we do not include any confining potential in 
our hamiltonian. Nevertheless, based on previous results 
obtained for the Hubbard model within the LDA or R- 
DMFT treatment including the trapping potential^ 
we expect our results to be reliable for the bulk of the 
optical lattice. The confinig potential typically affects the 
sites in a thin outer shell of the lattice, with a thickness 
of only a few sites^"^ . 

We now briefly discuss the DMFT method used to 
study the AHM. This method, which is exact in the limit 
of infinite lattice coordination, relies on a mapping of 



The retarded Green's function Qqo- is related to the full 
local Green's function Go- through the following self con- 
sistency conditions, 

g^J-iiUn) = iUn + A^CT - tlGa{iUJn) (4) 

LOn are the Matsubara frequencies and ^„ are the chem- 
ical potentials for the two species. At the self consistent 
point, the impurity and the lattice local Green's func- 
tions coincide. Therefore, the Eqs. (3) and (4) can be 
simply solved by direct substitution. Note that unlike the 
HM case, here the disordered phase is characterized by 
two self-consistent equations because the AHM explicitly 
breaks orbital rotational invariance for all r 7^ 1. 

2. Checkerboard-like ordered solution 

Given that both the HM and the FKM realize phases 
with long range order, we expect the AHM to also have 



the present problem onto a quantum impurity problem 
subject to self-consistency conditions. Here we consider 
the AHM on a Bethe lattice characterized by the semi- 
circular density of states, 

P.{e)^^yi~{e/D,Y (2) 

with Dc = 2tcr . We set = _D = 1 as the unit of energy 
in what follows. The advantage of the Bethe lattice is 
that it provides numerical simplifications and also is ex- 
pected to yield accurate estimates for a 3D simple-cubic 
lattice^, adopting D„ = 6t„. 

We solve the DMFT equations using two standard nu- 
merical methods: the Hirsch-Fye quantum Monte Carlo 
algorithm (HF-QMC), which is exact in the statisti- 
cal sense; and the exact diagonalization method (ED), 
which is in principle exact and whose discretization er- 
rors can be systematic reduced. Their implementation 
is described in detail in Ref. 2. The use of these two 
different methods further provides a crosscheck of our 
numerical results. 

To obtain a complete phase diagram of the AHM, we 
need to consider both ordered and disordered solutions 
to the dynamical mean field equations. We discuss the 
two cases below. 

1. Disordered solution 

In the disordered phase, the system can be mapped 
onto a quantum impurity problem described by the imag- 
inary time effective action: 



(3) 
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such phases. Based on the nature of the known solu- 
tions of the aforementioned models, here, we assume a 
checkerboard-like order for the AHM, as well. This im- 
plies the existence of two sublattices A and B such that 
each site of sublattice A is connected only to the sites of 
sublattice B and viceversa. In the AHM, the existence 
of two inequivalent sites necessitates a mapping onto two 
impurity models similar to the one described in (3) and 
the respective Green's functions obey the following self 
consistency conditions, 

GoaBi^^ri) = ioJn + l^-o- - 4G<xa(«W„) (5) 

In the limits r = 0, 1, this set of four equations reduces 
to two. For the case of general asymmetry (0 < r < 1), 
one needs to deal with the four different self-consistency 
conditions, which make the problem of the AHM con- 
siderably harder. Technically, we have to simultaneously 




solve two different quantum impurity problems, with two 
different effective baths each. Solutions with incommen- 
surate order which might exist in this model, especially 
away from half-filling^^'^*'^^ are not considered in this pa- 
per. From the structure of the mean-field equations Eqs. 
5, we find that for checkerboard-like order at half-filling, 
the Green's functions obey the following constraints: 

ReGAcr(it^n) = -ReGBCT(iwn) 
ImGAcr(it^n) = ImGBff(iwn) (6) 

We term this symmetry broken state as one with orbital 
order (00), and, in the symmetric case for r=l it coin- 
cides with the familiar Neel antiferromagnetic state. 

From Eqs. (5) and (6), we see that unlike the simpler 
case of the HM (r = 1) , here Gai and Gbh can be differ- 
ent for generic values of r^l. This implies the possibility 
of a charge density wave (CDW) order, as wc shall show 
below. 
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FIG. 1. Sketch of different ordered solutions of the AHM. 
Black (lower) dots, represent particles in the heavy band, 
while red (upper) dots represent those in the light-one. Big- 
ger dots represent higher occupation. In the SU(2) symmetric 
point r — > 1 (right-most pannels), only orbital order is possi- 
ble, and orbital polarization increases with U. As r decreases, 
CDW {ua 7^ ns) coexists with orbital order (lower left and 
central panel). CDW formation, nevertheless, decreases with 
increasing interactions (top panels), since in the Mott insulat- 
ing state, double occupied sites are energetically unfavorable. 



B. Long-range order phase 

Restricting ourselves to the half filled case with equal 
numbers of light and heavy particles, we find that in the 
orbitally ordered phase, the average occupations of dif- 
ferent orbitals on the A and B sublattices satisfy the 
constraint uai +nBi = UAh +nBh = 1- This implies that 

riAi - riAh = -{riBi - riBh) (7) 

Additionally, the sublattice occupation Uai + Uah 

{a = A, B) satisfies the relation 

- (1 - n^) = 1 - ns (8) 

Eq. (7) indicates that the orbital polarization in sublat- 
tice A is the opposite to that of B. The order parameter 
can now be defined as mo= {uai — Uah) with, a being 
either ^ or _B. In the HM limit, this orbital order trans- 
lates into a spin density wave with Neel antiferromagnetic 
order, and TOq is the familiar staggered magnetization. 

Note that Eq. (8) indicates that orbital ordering coex- 
ists with a charge density wave formation in which one 
sublattice has increased its occupation with respect to 
the other. The CDW order parameter can be defined as 
fnc= {nA — nB)/'2- We see that we have CDW order for 
all r^l. In the limit of the HM, both I and h orbitals 
arc equivalent, SU(2) symmetry is recovered and CDW 
order is lost, since uai = riBh and orbital order translates 
to Nccl antiferromagnetic order. An 00 with a CDW, 
in contrast, is realized in the FKM limit"*. On heating, 
these ordered phases survive upto a certain critical tem- 
perature To, which coincides with the Neel temperature 
in the HM, and with the Tcdw in the FKM case. To 
help visualize the variety of states, we sketch in Fig. 1, 
the different phases as a function of the interactions and 
the hopping asymmetry. 



III. RESULTS 



We now solve the relevant impurity problem with the 
associated self-consistency condition in both the ordered 
and the disordered phases, using exact diagonalization 
and quantum Monte Carlo methods^. The use of both 
methods allows us to crosscheck the results and also ob- 
tain the finite temperature phase diagram in the whole 
relevant parameter range of the model. In Fig. 2, we 
plot the finite-T Matsubara local Green's functions for 
different values of the parameters U and r =th/ti for 
T/D^D/ j3=Qmb. Although detailed features like the 
densities of states are technically difficult to extract from 
the Matsubara Green's functions, some basic properties 
can be readily obtained. For instance, the extrapolated 
value of G{iuJn — > 0) yields the density of states at the 
Fermi energy, which is expected to be finite in a metal, 
while it goes to zero in an insulator at low tempera- 
tures. This qualitative feature can be observed in the 
solutions for the disordered (ordered) phase shown in the 
left (right) panels of Fig. 2 for various values of r and 
the interaction U. This feature allows us to pinpoint the 
precise location of the metal-insulator transition in the 
model. Moreover, within DMFT, if the metallic solution 
is a normal Fermi liquid, it can be shown that the den- 
sity of states at the Fermi energy at T = for any U, 
remains pinned at the U — value^. Consequently, if the 
extrapolated value of the Green's function does not go to 
zero (i.e., is not an insulator), but violates the pinning 
condition, then the self-energy would be non-vanishing 
at a;„ — s- 0, hence the states have a finite lifetime, and 
the solution corresponds to a non-Fermi liquid metallic 
state. All of these features, coupled with the estimation 
of various order parameters, help us obtain an interesting 
and rich finite temperature phase diagram for the AHM. 
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Ordered 




FIG. 2. ImG(ia;„) for T/D = 0.025. The left (right) col- 
umn corresponds to the disordered (ordered) solution for 
U = D,2D,3.5D, for the first, second and third rows, re- 
spectively. Empty (filled) symbols correspond to the light 
(heavy)-orbital, while circles, squares and diamonds corre- 
spond to r = 0.1,0.4, 1, respectively. For the ordered state, 
note that as shown in Eq. 6, ImG(ia;„) is independent of the 
sublattice. 



A. Disordered phase 



Phase 



The disordered phase of tlie AHM, lias been discussed 
in detail in recent publications (Refs. 16 and 21). Here, 
for completeness and to facilitate the ensuing discussion, 
we briefly summarize the main features of the disordered 
phase diagram. It is schematically shown in Fig. 3 for a 
generic value of the mass asymmetry parameter r. This 
phase diagram would be appropriate for a fully frustrated 
system. 

At low T, there exists a metal-insulator transition from 
a metal to a Mott insulator as one increases the repulsion 
U. The metallic phase is characterized by a Fermi liq- 
uid for all r 7^ and generically has two kinds of quasi- 
particles with different masses. The light particles are 
more strongly renormalized, since the screening poten- 
tial produced by the heavy ones is less efficient due to its 
lower mobility. Despite the different mass renormaliza- 
tion, both quasiparticles disappear (their mass renormal- 
ization diverges) at a single critical U = Uc2, signaling 
the metal-insulator transition. On the other hand, if one 
starts at high values of the intereaction U, well in the 
Mott state, and progresively decreases the strength, the 
insulating gap decreases, and eventually closes. At that 
point there is an insulator to metal transition, which is 
denoted as Uci- In between these two critical values, one 
finds a coexistence regime, with two solutions, one in- 
sulating and one metallic, as depicted in Fig. 3. The 
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FIG. 3. Schematic T — U phase diagram for a generic value of 
the parameter r^^ . At low T there are three distinct regimes: 
a Fermi-liquid state with two different quasiparticles; a co- 
existence zone, and a Mott insulating phase. In the metallic 
regime, above Tcoh there is a crossover to a Falicov-Kimball- 
like non-Fermi liquid state. Increasing interactions from this 
phase, the system crosses over to a bad insulating state (a 
finite T Mott insulator). 



coexistence region extends in temperature from T^O up 
to finite temperatures. Both solutions are thermodynam- 
ically stable within the disordered phase, thus the metal- 
insulator transition is first order at finite T and culmi- 
nates in a second order critical point, similarly to the 
HM case^ . The critical values of the interactions and the 
T- value of the critical end point of this regime, i.e., the 
area of coexistence, decrease with decreasing r. 

At lower U, the metallic Fermi liquid state survives 
upto a coherence temperature TcohiU,r). Above Tcoh, 
the system crosses over to an orbital selective metallic 
phase, which closely resembles the corresponding non- 
Fermi liquid state of the FKM model. This intermediate 
temperature regime, can be thought of as if the hopping 
of the heavy particles were effectively renormalized to 
zero, hence becoming localized, while the light ones re- 
main mobile. Orbital selective crossovers between quali- 
tatively different metallic phases, though commonly ob- 
served in real materials, such as high-Tc pnictides, heavy 
fcrmions, etc., are often not well understood. In this con- 
text, our model may be viewed as a minimal model for a 
non-Fermi liquid state emerging from an orbital-selective 
mechanism. 



2. Momentum dependent density of states 

Since the AHM is more generic than the Hubbard 
model, it may be easier to realize in cold fermionic mix- 
tures. Moreover, from the previous discussion it should 
be clear that the AHM presents the tantalizing possi- 
bility of realizing non-Fermi liquid states at relatively 
high temperatures. These types of poor metallic states 
arc central to many open problems in condensed matter 




FIG. 4. Density distributions of the light and heavy orbitals 
for r = 0.2 and U/D = 0.8. As T decreases, n{e) tends to the 
usual Fermi-liquid distributions where there is a discontinuity 
at the Fermi energy (black arrows show the tendence of n(e) 
as T reduces). The vertical black line at e = shows the 
expected jump at T = 0. 
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FIG. 5. Slope of the light particles density distribution as 
a function of [T / D)'^ , for U/D = 0.8. For T > Tos, the 
metallic behavior of the AHM, coincides with the non-Fermi 
liquid metallic state of the FKM. For this value of U , Tos(r = 
0.4) ~ 0.2D and Tos(r = 0.2) ~ O.ID. 



B. Ordered phase 



physics. This immediately raises the question of what 
would be a measurable physical quantity, which is sen- 
sitive to the orbital selective crossover to the non-Fermi 
liquid state. 

One practical probe would be time-of-flight 
experiments^, where the momentum resolved den- 
sity distribution n{k) is measured. In a conventional 
Fermi-liquid at zero temperature, n{k) has a jump of 
size Z < 1 at the Fermi momentum kp, with Z being 
the quasiparticle residue. Replacing the momentum by 
the single particle energy e, the resulting quantity n(e), 
which is more easily obtained within DMFT, would also 
exhibit the same jump at the Fermi energy e ~ ep. At 
nonzero T, the discontinuity in n is smeared within an 
energy scale given by T (sec Fig. 4). In a non- Fermi 
liquid state, Z = 0, and there is no discontinuity even at 
T = 0. 

Since the crossover of the AHM is at finite T, we do 
not expect to see sharp changes in the density distribu- 
tion. Nonetheless, going from the HM towards the FKM 
(reducing th/U and keeping all other variables constant), 
n(e) evolves to n^^'^(e), for r < Tcrit- However, we find 
that the slope of n;(e) at the Fermi energy, as a function 
of is a better magnitude to pinpoint this crossover. 
We plot the slope in Fig. 5 and we clearly see that above 
a certain temperature Tqs, which depends on r and U, 

AHM cj FKM 

^^^^gj — approaches ^^^^ — . Therefore, this signals that 
the light orbital of the AHM evolves to the non-Fermi 
liquid state of the FKM once orbital-selective crossover 
has been achieved. 



In the non-frustrated case, we expect the system to 
order at low temperatures. Solving the equations dis- 
cussed in section II A 2, we obtain the critical tempera- 
ture To{r, [/), below which the orbital order sets in, and 
the behavior of the order parameters TOq and mc, defined 
before. 



1. Ordering Temperature To 

We find that analogous to the HM and FKM, at T = 
and for all [/ 7^ 0, the system is an insulator with long- 
range order which vanishes above a critical temperature 
To(U,r). The results are shown in Fig. 6 for a generic 
value of r=0.4. For comparison, in the figure we also 
show the main features of the phase diagram discussed 
before, in Fig. 3. We observe that the metal-insulator 
transition, described in section HI A is shadowed by 
the 00 phase, which is thermodynamically more stable. 
This behaviour is similar to the one realized in the Hub- 
bard modeP. It is interesting to observe that both the 
coexistence region, the orbital-selective thermal crossover 
and the maximum in the orbital ordering temperature To, 
all occur at the same parameter region, where the inter- 
action U is comparable to the model bandwidth 2D and 
the correlation effects are strongest. 

Within DMFT, the simplest procedure to compute the 
critical temperature To is to start at low T, where the 
converged DMFT solution is found to have orbital long- 
range order (i.e., ruo^ 0, see eq. 7) and then slowly in- 
crease the temperature until the order vanishes at some 
T = To{r,U). We note that these transitions though 
continuous and second order, are rather sharp as can be 
seen in Fig. 7, that shows the order parameter mo(T), 
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in the case of the HM^. The different phases of the disordered 
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FIG. 7. Orbital order parameter mo for different values of the 
asymmetry r and U — D. 



for U ~ D and several values of the hopping ratio r. In 
Fig. 8, we plot the ordering temperature as a function of 
both, the interaction U and hopping asymmetry r. 

At weak coupling i.e., small U and for all ?- ^ 0, wc 
see an exponential behaviour , To oc exp(— c/J7), where c 
is some constant, similar to the Hartree mean field result 
for the Hubbard modcl^®. In the FK limit, we expect^ 
To oc (U/D)^ \n{D/U). The qualitative dependence of T^ 
on r can be understood by considering the competition 
between kinetic and Coulomb energy. In the disordered 
metallic region, when [/ = 0, To goes to zero, and there 
is no order. In the HM, as U increases, the ratio between 
potential and kinetic energy (Ep/Ek) increases, and so 
docs To.^ Since a decreasing r decreases Ek, the ratio 
EpI Ek and hence To increase as r decreases in accordance 
with Fig. 8. 

To understand the non-monotonic behaviour of To seen 
in the Mott insulating side of the phase diagram, it is use- 
ful to consider the large-?/ limit, where the AHM can be 
mapped^® onto a spin-1/2 anisotropic Heiscnberg (XXZ) 




FIG. 8. Critical temperature To as a function of interaction 
U (upper panel) and hopping ratio r (bottom panel). The 
scaling To oc (1 + r^) at strong coupling, is shown in the inset 
of the upper panel. 



model described by the effective hamiltonian, 



a spin — J / ^ [ 



-7K^|-1)] (9) 



with exchange constants J =tith/U and 7 = ^j^- 
This model interpolates between the isotropic Heisenberg 
model at r = 1 and the Ising model at r = 0. A mean 
field study of this spin model shows us that the ordering 
temperature To in the large U limit is proportional to 
J-f (X {1 + r^) consistent with the DMFT results, which 
show that To decreases with decreasing r in the large U 
limit (see inset of Fig. 8). This feature also persists in 
the intermediate U regime as shown in the lower panel of 
Fig. 8. Therefore, from the previous discussion, the non- 
monotonic behaviour of Tq seen in Fig. 8 is a consequence 
of the increasing dependence of the ordering T with U at 
weak coupling, on one hand; and of the decreasing de- 
pendence at strong coupling due to the super-exchange 
interaction, on the other. 



2. Orbital and charge order 

We now discuss the orbital and CDW ordering seen in 
the ordered phase for T < To- In Fig. 9, we plot the or- 
der parameter TOo (defined in sec. II B) as a function of 
the asymmetry r for different values of the interaction, at 
a fixed finite low T. The variation of mo with r is rather 
complex and depends on the value of U. For large enough 
U, we observe a sudden increase of the orbital order pa- 
rameter with increasing r. This is because the energy of 
the 00 state at large U, where we have a Mott insulator, 
is given by the super-exchange, which increases with r, 
as discussed above. At lower values of [/, in contrast. 
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the system is in a more itinerant state, a Slater like 00. 
Thus, increasing r the system is driven by kinetic energy 
towards a metallic state, lowering the magnitude of the 
order parameter. From a practical standpoint, these re- 
sults suggest that in experiments, for any given value of 
the interactions U, the phase transition to the 00 phase 
can be attained by tuning either T or r. This second 
route to achieve the long-range ordered phase can be of 
particular interest in cold atom systems where it is often 
very difficult to reach the low temperatures associated 
with long range ordering. 

Another interesting aspect of the ordered phase, is that 
orbital order coexists with a genuine charge density wave 
(CDW) state. This is illustrated in the bottom panel of 
Fig. 9 where the CDW order parameter rric (defined in 
sec. II B) is plotted as a function of the hopping asymme- 
try for different values of the interaction. We observe that 
the total number of particles per site is different depend- 
ing on whether the site belongs to the A or B sublattice, 
as qualitatively depicted in Fig. 1. The behavior of nic 
is non-monotonic, large at low interactions and small r, 
and smaller as one approaches the Hubbard limit, where 
it eventually vanishes. This can be understood from the 
fact that large U is detrimental for charge inhomogene- 
ity, and penalizes the formation of a CDW. Alternatively, 
an increase of r also favors a more homogeneous state, 
reducing the order parameter mc- The situation at low 
U is qualitatively different and can be understood using 
the same arguments presented earlier for the behavior of 
mo- 

Further insight is gained by computing the contribu- 
tion of each orbital to the formation of the CDW. In 
Fig. 10 we plot the mean number of particles in each 
of the two orbitals on sites of the A and B sublat- 
tices. We find that the occupation of the heavy orbital 
has a stronger staggered character than the light one 
i\nAh " nBh\ > \nAi — n-BiD- This can be understood 
by a continuity argument: in the FKM limit (r=0) and 
low U, the heavy orbital is localized while the light one 
remains itinerant. For higher values of U, where one is 
within the ordered phase for all r, we see that this pro- 
nounced staggering persists, though becomes weaker as r 
increases and eventually vanishes at the Hubbard point 
(r=l). 

All of these results are consistent with the behaviour of 
the double occupancy, which can actually be experimen- 
taly measured in cold atom systems^^. In Fig. 11, we 
plot the double occupancy d — {nhfii) in both sublattices 
in the ordered phase and also its value in the disordered 
phase. The fact that true CDW order accompanies the 
orbital order immediately implies that d can show differ- 
ent behaviours on each of the sublattices. This is indeed 
the case as shown in Fig. 11, where dA and ds are dif- 
ferent within the ordered phase (T < To), when r ^ 1. 

In the symmetric case, r = 1, we observe that, for mod- 
erate to strong interactions, the double occupancy of the 
ordered phase is higher than in the disordered one. This 
is due to the fact that in a state with staggered order. 
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FIG. 9. Orbital order parameter (upper panel) and charge 
density wave order parameter (bottom panel) as a function 
of r for different values of U and T /D = 0.05. Inset: To for 
different values of r as a function of U . The liorizontal dashed 
line corresponds to T/D = 0.05, which is the temperature 
used for computing mo and rric. 
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FIG. 10. Left panel: Number of light and heavy particles on 
the A and B sites as a function of r for different values of 
U and T/D — 0.01. Right panel: ua and ub- Notice that 
7^A 7^ nB for U—D, while for U—0.2D the CDW disappears 
for r > O.ID. 



hopping to a neighboring site is less penalized by Pauli 
blocking with respect to the situation in a disordered 
state^'^. On the contrary, for small to moderate coupling 
the behavior is the inverse, since in the disordered metal- 
lic state the particles delocalize on both orbitals, increas- 
ing the probability of double occupation. The ordered 
state remains an insulator even at small U, thus favor- 
ing orbital localization. At low enough U, the ordered 
state eventually becomes unstable at the finite T of our 
numerical calculations, ie, when T > To{U). 

For r 1, we observe the effects of the onset of the 
CDW. The behavior of the double occupancy becomes 
different on each sublattice, as previously discussed. At 
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FIG. 11. Double occupancy (Ia and ds (solid lines) for dif- 
ferent values of U and r at T/D = 0.01. For comparison, 
the dashed lines represent the equivalent curves for the dis- 
ordered state. The sharp behavior at low f/ is a signature of 
the disorder-order transition at finite T. 



low couplings, we find that the difference has a sudden 
increase, which reflects the steep rise of the orbital order- 
parameter Too as the system enters the ordered phase. 
This shows that an abrupt change in the sublattice dou- 
ble occupation is a good probe of CDW order in a sys- 
tem. This sharp behaviour reduces with increasing r as 
expected because the CDW order decreases with increas- 
ing r, fully disappearing in the limit of r = 1. We also see 
that, for any value of U, the double occupancy is higher 
in the HM and decreases with decreasing r which is con- 
sistent with the fact that as r decreases it gets harder 
for the heavier particle to hop from one lattice site to 
another. These results confirm that asymmetric hopping 
amplitudes favor a CDW, which can be experimentally 
detected. Alternatively, this also suggests that the den- 
sity distribution can be manipulated by tuning the hop- 
ping ratio r. 



3. Thermodynamics 

One of the challenging aspects of cold atomic systems 
in optical lattices is the determination of their temper- 
ature. To this end, it is important to know the ther- 
modynamic properties of the system. With this issue in 




FIG. 12. Specific heat and entropy per site as a function of 
T for r = 0.4, in the disordered case for U = D (metal) and 
U = 2D (insulator). For reference, Tcoh is indicated by the 
vertical dotted line. It corresponds to the destruction of the 
Kondo effect, the first peak in the Cy and the ln(2) plateau 
in S/N. 



mind, we focus on the computation of the fundamental 
thermodynamic quantities: specific heat and entropy of 
the AHM in this section. We set the Boltzmann constant 

In our previous paper^^, we discussed the specific heat 
Cv and the entropy per site (or equivalently per particle 
if working at ntot = 1) S/N=s in the disordered phase 
of the AHM. The main features can be summarized as 
follows (see Fig. 12 where we plot these quantities at a 
finite generic value of r): In the metallic phase at low 
U, the Cv and s are linear in T, for all values of the 
asymmetry parameter r ^ and T < Tcoh- This be- 
havior is a signature of the Fermi liquid groundstate. At 
T « Tcoh, the Cy develops a peak and the entropy has a 
small plateau at s = ln(2), that corresponds to the two 
orbital degrees of freedom. In the Kondo impurity prob- 
lem, this scale indicates the destruction of the Kondo 
singlet state [{\lh) - \hl)) (|0® I''-))] • Increasing T 
further, s increases and eventually saturates at s = ln(4) 
at high T. This value corresponds to the four different 
states available per site, ie |0), \h), \l), \lh). On the other 
hand, Cy develops a second peak due to the excitation 
of states of the incoherent Hubbard bands. 

The behavior of the entropy in the Mott state of the 
AHM, shown in Fig. 12, resembles that of the metallic 
state of the standard FKM^^. However, the similarity is 
superficial, and is due to different physical reasons. In the 
latter, the heavy particles have zero hopping amplitude 
and are localized, thus there is an equal probability of 
finding a heavy particle or not on any given site. There- 
fore, there are two occupation states associated with any 
site and this contributes a factor of ln(2) to the entropy, 
even as T ^> 0. As T increases, the entropy increases 
monotonically to the asymptotic value ln(4). However, 
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FIG. 13. Specific heat and internal energy as function of T 
for U = D, and different values of r in the ordered phase. The 
first peak in Cv reveals a formation of an orbitally ordered 
phase, while the second, the localization of the particles. The 
sharp reduction of the internal energy when decreasing T, 
reveals the formation of an insulated-orbitally ordered state, 
at To. For U — D, To is weakly dependent on r, as shown in 
Fig. 8. 



in the Mott insulating state of the AHM, both particle 
species are localized and, in contrast to the FKM case 
that we just described, there are two different ways of 
singly occupying a site or |/)), hence the ln(2) in the 
entropy. As before, at high T, density fluctuations arc 
thermally activated, hence the entropy and s also goes 
continuously to ln(4) . Regarding the behavior of the Cy 
since the insulating Mott state of the AHM is character- 
ized by the opening of a correlation gap, the specific heat 
displays activated behavior at low T. 

The behavior of the specific heat in the ordered phase 
shows clear signatures of a phase transition. The results 
arc shown in Figs. 13 and 14 for characteristic values 
of the parameters r and U. At very low T, the Cy has 
activated behavior (the internal energy is almost con- 
stant), due to the opening of the insulating gap in the 
orbitally ordered state. Increasing T, the system even- 
tually reaches the transition temperature T = To where 
the order is lost. The internal energy plotted in Fig. 13 
shows a sudden upturn as the state becomes unstable. 
At the critical point there is a kink which translates to a 
divergence in the specific heat. 

In Fig. 15 we plot the entropy in the ordered state 
for r = 0.4 and several values of U. The data shows ac- 
tivated behavior at very low T and a rapid increase of 
the thermal excitations as the temperature approaches 
To and the ordered moment melts (see Fig. 7). For T 
above To, in the disordered phase, the entropy is iden- 
tical to that of Fig. 12. The non-monotonic behavior 
of the position of the maxima, is simply a consequence 
of the peculiar behaviour of To with U that was already 
discussed Sec. IIIBl. 
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FIG. 14. Specific heat and internal energy as function of T 
for U — 2D, and different values of r in the ordered phase. 
Note that for U = 2D, To increases with r. 
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FIG. 15. Specific heat and entropy as a function of T for r = 
0.4 and U/D = 1, 1.3, 2,3, in the ordered phase. Above To((7), 
where the Cv has a kink, the entropy coincides with that of the 
disordered state. For low enough T, the activated behavior 
is non-monotonic in U, revealing the different natures of the 
ordered state with varying U. 



C. Proposal for cooling mechanism: Entropic 
chromatography 

Despite the fact that experiments on cold atoms in 
optical lattices are done at ultra low temperatures, it is 
necessary to achieve even lower temperatures, to access 
new interesting strongly correlated quantum phases. 

One of the theoretical schemes proposed for cooling, is 
based on the idea of embedding a fermionic gas in an op- 
tical lattice within a bosonic bath, and then 'squeezing' 
out the entropy of the fermions. This may be achieved by 
first shrinking their trap, thus compressing the fermionic 
gas towards a band insulator state. This state has no 
free degrees of freedom left, hence, no entropy (with the 
exception of the states at the edges of the trap). The 
lost fermionic entropy is thus transferred to the bosonic 
bath, and the bosons are then released. The remain- 
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ing fcrmions are in a low entropy state, therefore at a 
much lower temperature'^''. Another proposal ^, sug- 
gested modifications to the shape of the harmonic trap 
so as to create a deep dimple in the center where one 
would form a band insulator. This separates the system 
into an entropy poor state at the center, surrounded by 
entropy rich regions. The entropy rich particles could 
subsequently be removed from the system by partially 
opening the trap and the remnant band insulator could 
then be adiabatically relaxed to the relevant experimen- 
tal regime. 

Motivated by the present study of the AHM, here, we 
propose an alternative scheme to cool the system based 
on the principal feature of the AHM i.e., the hopping 
asymmetry. In Fig. 16, we plot the entropy per site s, in 
the non-interacting limit (C/ = 0) limit of the AHM at 
half filling. We see that most of the entropy is carried 
by the heavier particles for T < 0.4D. The difference 
in entropy between the two species is more pronounced 
at very low temperatures. Though we have not included 
the effect of the trapping potential in our calculation, we 
do not expect it to significatly change the nature and 
validity of the cooling mechanism that we shall describe. 

The ability to concentrate most of the entropy in the 
heavy species by manipulating the hopping amplitudes of 
one species, immediately suggests the possibility of an en- 
tropic chromatography method. One can repeatedly con- 
centrate entropy in one of the species and use the other 
as a thermal bath, which is then partially evaporated, 
liberating entropy (ie heat) at each successive step. This 
leads to a big reduction of entropy of the remanent sys- 
tem, hence reducing its temperature. Our cooling scheme 
is based on this cntropic chromatography, supplemented 
with a series of system parameters modifications, to en- 
sure that the physical system stays within the required 
filling regime (ie one particle per site on average). This 
scheme is similar to those mentioned above, in that we 
separate the system into entropy rich and entropy poor 
sectors, with the main difference being that here the cn- 
tropic separation is not spatial. The advantage of this 
procedure is that it does not require an additional Bose 
condensate to act as a reservoir or additional lasers to 
modify the trapping potentials in a precise manner. Our 
scheme essentially relics on the possibility of having two 
different fermionic atomic species within their respective 
atomic traps, which can be independently adjusted. We 
describe our cooling procedure in detail below. 

We consider a system where the light and heavy par- 
ticles have their respective trapping potentials and their 
interaction is described by the AHM. We limit our anal- 
ysis to the disordered phase of the AHM. However, the 
extension of the idea to the long-range ordered phase can 
be performed along similar lines. In this section, for the 
sake of clarity in the description of the procedure, we de- 
note the two atomic species as particles of type 1, and 
type 2. We assume that all the system parameters (hop- 
ping, U, trapping potential) can be tuned adiabatically 
and that the system thermalizes to an equilibrium state. 
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FIG. 16. Entropy per lattice's site in the non-interacting sys- 
tem for r — 0.2. 



The cooling scheme to be described below is useful and 
applicable provided the initial temperature and entropy 
of the system is such that in the non-interacting limit, the 
difference in the entropies of the light and heavy particles 
is substantial. 

Wc summarize the cooling scheme in the following 
steps : 

i) Start by loading Ni ^ N2 ^ N particles of each 
fermionic species into the optical lattice at finite U 
and hoppings 



t2 



t 

at 



(10) 



with an asymmetry a <^ 1. The parameter a ac- 
tually plays the same role as the ratio r, however, 
we use this convention to emphasize that a is varied 
during the cooling procedure. The entropy of this 
initial state is denoted as S).Ji. The trap and the 
optical lattices are adjusted such that the system is 
half-filled. 

ii) The interaction is now adiabatically tuned toU = 0. 
Since this is an isentropic process, the entropy of the 
resulting system is the same as that of the initial 
state. However, as shown in Fig. 16, in this non- 
interacting limit, the entropy of the system is now 
mostly carried by the heavy particles 2. 

For non-interacting fcrmions (or quasiparticlcs in a 
Fermi liquid), at low T the entropy is proportional 
to T^"^: S = K^kprnT /iTi where the mass of the 
particles (or the effective mass of the quasiparticles 
in the Fermi-liquid regime) m is proportional to the 
inverse of the hopping amplitude t and kb is the 
Boltzmann's constant. This implies that at this step 
of the cooling process, the entropy of the system can 
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be expressed as: 



c.{») 
'-'tot 



Ns 



(1) 



Ns. 



(1) 



at 



1 

1 + - 

a 



(11) 



where s^^^ is the entropy per particle of species j=l, 
2 at this step, and T^^-* is the new temperature re- 
sulting from the adiabatic variation of U. k is some 
proportionality constant which is the same for both 
species (cf. Fig. 16). 

iii) Open the trap of particles 2, and free half of them. 
The released entropy (5*1^^) is: 



:-(!) 



2at 



(12) 



In order to compare the temperature of this final 
state with the initial one, we turn off U. This fi- 
nal state has equal numbers of particles N/2 of each 
species and the same hopping ratio as in step i and 
an entropy 



q(f) 

'-'tot 



c(3) 
'-'tot 



KNT^f^ 
2t 



(1/a + l) 



(16) 



The isentropic processes also ensure that the temper- 
ature of this state T^^'^ = T*^^^ , which follows from 
Eqs. (14), (15) and (16). Then, by matching the final 
entropy of Eq. 16 with 



:-(/) 



S, 



(1) 



we find that 



TU) _ 2a + 1 
TW ~ a + 2 



(17) 



(18) 



iv) Turn on U adiabatically so the particles interact and 
thermalize (we assume here that the relaxation time 
is fast). The system now has an unequal number of 



particles of both species and a total entropy S. 



(2) 



'-'tot 



(1) 



v) To eliminate the excess number of particles of species 
1 and to reduce the entropy further, we adiabatically 
invert the ratio ^2/^11 h^ order to have. 
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at 
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vi) The interaction is again adiabatically tuned to J7 = 

(2) 

0. We now have a state with total entropy S^^l such 
that 



j(2) _ 
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where sp'' and ' denote the entropies per particle 
of the two species and T^^^ the resulting temperature 
at this step. 

vii) We open the trap for particles 1, and eliminate N/2 
particles of type 1. The freed entropy at this step is 



,(2) 



^(2) 



KNT/f> 
2at 



and the total residual entropy is 
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viii) To obtain the final state with N /2 particles of each 
species and to thermalize the system so as to have a 
well defined temperature, we, once again, invert the 
ratio ^2/^1, at finite U . The traps of both species 
are narrowed, adapting the number of sites ensuring 
that we have a half-filled system. 



Note that the initial temperature of the mixture T*^*' is 
not know, the first temperature that could be defined was 
T(^) at step (n) with Eq. 11. From the last expression 
and Eqs. 11 and 16 it also follows that the ratio of the 
total final and initial entropies is s[l] / sf^t = T^f) /2T'^^\ 
thus that the entropy per particle in the f7 = limit can 
be reduced by the same factor r(/)/T(i). Wc also see 
from (18), that no reduction of entropy or temperature 
is achieved in the Hubbard limit a 1 and that the 
cooling is entirely a consequence of the anisotropy in the 
hoppings. The maximal reduction T^-^^ = T^^^ /2 is ob- 
tained in the FK limit, a — ?► 0. Note that depending on 
the number of particles loaded into the optical lattice and 
the final system required, the cooling procedure can be 
repeated a finite number of times, resulting in a drastic 
reduction of the temperature. 

However, we should mention that in cold atom exper- 
iments, currently accessible entropies per particle are of 
the order of s ~ 1.5. This value is well outside the lin- 
ear regime discussed here (cf. Fig. 16) and the analyti- 
cal result obtained here, though helpful to illustrate the 
mechanism, is quantitatively invalid. Nevertheless, we 
emphasize that the cooling procedure i-viii, should still 
qualitatively work within a wider window of T, since the 
large difference in entropy persists well beyond the linear 
regime (cf. Fig. 16). 

Within our approximation, we now attempt to esti- 
mate the efficacy of our cooling procedure for the inter- 
acting system. In Fig. 17, wc plot the entropy of the 
interacting ([/ =D) and the non-interacting system, for 
^2/^1=0. 2 at half-filling. In this example with q;=0.2, we 
see that starting from U=D and entropy per site (at half 
filling is identical to per particle) s=l, the initial temper- 
ature is about T^^^/D ~ 0.28. Then, first by decreasing 
U to zero, we see that it corresponds to a temperature 
r(i' ~ 0.1751). Following the steps i-vm, we can use 
eq. 18 to estimate the temperature of the new state. In 
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FIG. 17. Entropy as a function of T with a hopping asym- 
metry t2/ti = 0.2, for [7 = and U = D. T^^'l {T^^{) corre- 
sponds to the temperature of the non-interacting (interacting) 
state, with S/N = 1, while T'-'' {Tj;J^t) correspond to that of 
the final state (after doing the cooling procedure, and taking 
a = 0.2). 



this example, the temperature of the non-interacting sys- 
tem would be T'-f'^ ~ 0.1 ID, which according to Fig. 17, 
corresponds to an entropy of s ~ 0.83. Then, adiabat- 
ically turning on the interactions to U=D, it would be 
possible to get an interacting state that corresponds to 
a temperature of T>^1 / D ~ 0.16, which is about a 40% 
lower than the initial temperature However, we reiterate 
that this number is only an estimate, since we have used 
Eq. 18 which is accurate only in the linear regime. In 
any case, the percentage could be further increased by 
increasing the hopping asymmetry, i.e., reducing a. The 
cooling procedure can be repeated a number of times to 
access even lower entropies and temperatures. It would 
also be interesting to see how this scheme works if we take 
into account the trapping potential and also consider or- 
dered initial states. This can be done numerically and is 
left for future work. 



Regarding the actual implementation of the cooling 
procedure in experimental setups, we would like to point 
out that our proposal has the simple goal of showing that 
the hopping asymmetry can potentially be exploited to 
cool the system once it has been loaded onto an optical 
lattice. Though the hopping and interaction parameters 
can in principle be tuned over a wide range, a main exper- 
imental concern is how to remove the excess of entropy 
from the system in a controlled manner. This is a very 
challenging issue and is the subject of ongoing work. An- 
other important aspect is the equilibration times of the 
system. These are topics of great experimental relevance 
and remain beyond the scope of the present study. 



D. Cooling by tuning hopping and interactions 

The success and applicability of the cooling procedure 
discussed above depends crucially on the changes in tem- 
perature induced by the adiabatic variation of the model 
parameters in the AHM. More precisely, a key issue is 
whether adiabatic changes of U and will heat or cool 
the system'^'^''^''^. In this section, we analyze the tempera- 
ture changes in the AHM as the parameters U and r are 
varied. As in the previous section, we restrict our study 
to the disordered phase, but the generalization to a long- 
range ordered phase is straightforward. From Refs. 33 
and 34, wc know that any change of a parameter x in 
the Hamiltonian, will result in an entropy change which 
is linked to a temperature change as follows: 



1 5s 
c 5x 



1 5T 
T~5x 



(19) 



where c = T5S/5T is the specific heat (c is always possi- 
tive). 

Using Maxwell's relations, 



5s_ 

5U 



N,t^,T 



5d 
5T 



5s_ 

5t: 



1 5Ek^. 



t„ 



5T 



(20) 



(21) 



and (19), we obtain the following expressions for the 
change in temperature for fixed entropy, when the Hamil- 
tonian parameters t/, are varied. 



1 5T 
T 5U 



1 5d 
c Jf 



1 5T 
T SL 



S,N,U 



1 5Ekin,cr 



5T 



(22) 



(23) 



N,t„.,U 



The changes in temperature when the interaction U is 
varied is given by the temperature variation of the double 
occupancy, whereas, the change in temperature induced 
by the variation of i^r is dictated by the derivative of the 
kinetic energy of the particle species a with respect to 
the temperature. As explained in Ref. 34, Eq. (22) 
suggests that increasing the interactions will cool (heat) 
the system whenever the temperature derivative of the 
double occupancy is negative (positive). In a similar way, 
Eq. 23 shows that increasing the hopping amplitude will 
heat (cool) the system if the temperature derivate of the 
kinetic energy is positive (negative). 

We now analyze the consequences of tuning theses pa- 
rameters in the particular case of half-filling without the 
trapping potential. In Fig. 18 we plot the double oc- 
cupancy as a function of T for various values of the in- 
teraction U and r. As pointed out in Ref. 33 for the 
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FIG. 18. Double occupancy as a function of T, for four values 
of the interaction (one for each panel), and three different 
values of the ratio r. 



Hubbard model, at weak coupling, increasing JJ heats up 
the system if one is at high enough T because 5d/5T > 0. 
For T < T„i, where Tm is the temperature at which the 
double occupancy is the lowest, 6d/ST < , we have a 
Pomeranchuk like effect where the system cools with in- 
creasing U. This cooling effect is lost at strong couplings, 
where the system is insulating and the double occupancy 
does not vary at low T and always has positive slope. 
Consequently, there will be neither cooling nor heating 
of the system whenever T is small, and there will be heat- 
ing if T is big enough to fill the gap. We find that all 
of these features, first predicted for the HM, exist in the 
AHM case also (cf. top panels of Fig. 18). 

We now focus on temperature changes induced by the 
variation of the hoppings. In Fig. 19, the kinetic energy 
for each species is plotted as a function of T. For low 
U, i.e. in the metallic phase, the kinetic energy of both 
species increases monotonically with T for all r. From 
Eq. 23, this implies that there will always be heating 
of the system as one increases to- in the entire temper- 
ature range of interest. In the insulating regime, since 
the system is in a Mott state, the kinetic energy remains 
practically constant at low T for both species (especially, 
the heavier one). Despite the minor variations seen at 
higher T we expect that a change of hopping parameters 
at large U leaves the temperature of the system effectively 
unchanged. These results bode well for the cooling proce- 
dure outlined earlier. It would be very interesting to use 
the results of this section to calculate the exact change in 
temperature resulting from finite variation of the Hamil- 
tonian parameters during the cooling cycle. This is left 
for future work. 



IV. CONCLUSIONS 

In the present work we have studied the AHM as a 
minimal model to explore multiorbital correlated physics 
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FIG. 19. Kinetic energy as a function of T for i-particles 
(left panels) and /i-particles (right panels). Upper panels: 
U/D = 0.5; lower panels: U / D = 2 (filled symbols) and 
U/D = 3.5 (empty symbols); black circles, red squares and 
blue diamonds corresponds to r = 0.1, 0.4, 1, respectevely. 



in cold-atoms. We obtained the phase diagram includ- 
ing both ordered and disordered phases. The disordered 
phase was found to exhibit an interesting orbital selective 
crossover at finite temperatures to the non-Fermi liquid 
phase that we first described in Ref. 21. We showed that 
the slope of the momentum dependent density distribu- 
tion at the Fermi energy may be a convenient observable 
to probe and characterize this orbital-selective state at 
finite T. In the ordered phase, which was the main fo- 
cus of the present study, we showed that the AHM dis- 
plays orbital order cocxisiting with a charge density wave 
whenever r ^ 1. We computed the behavior of experi- 
mentally accesible observables, such as double occupancy, 
which provided clear signatures of the onset of the or- 
dered state. We also found the regimes of hopping ratio 
and U where the ordering temperature is highest, that 
is, where it may be easier to access the ordered phase. 
We also discussed the relation between orbital order and 
the onset of charge density waves, which may provide a 
practical route for manipulating the density distribution 
inside optical lattices. 

Finally, we presented a new strategy to cool fermions 
in optical lattices using entropic chromatography exploit- 
ing the features of entropy concentration that follows 
from hopping asymmetry. In fact, currently achievable 
entropies per particle s ~ 1.5 are still too high to even ob- 
serve the intermediate temperature orbital selective state 
discussed in this paper, which is characterized by an en- 
tropy per particle s ~ 0.5 — 0.8. This cooling proce- 
dure might be a very efficient method to achieve the low 
temperatures required to access exotic correlated states 
in cold fermionic systems in optical lattices. An obvi- 
ous extension of our work would be to include the trap- 
ping potential in the AHM. Though this would result in 
quantitative changes to the the ordering temperatures, 
entropies per particle in the various phases, etc., we still 



15 



expect the phase diagram and the physics discussed in 
this paper to be vahd in the bulk of the system. We 
hope that this work may motivate further experimental 
developments and provide new bridges between the area 
of cold atoms and condensed matter physics. 

E.A.W. acknowledges support from the French pro- 
gram ANR-09-RPDOC-019-01. 



For the ordered phase, this expression can be generalized 
to 



Skin,. ^fJ2 GA<.(*w„)GB.(iw„)e'""°^ (A5) 



Appendix A: Analytical expressions of the 
observables 



For completeness, in this appendix we present the dif- 
ferent expressions used to compute the physical observ- 
ables from the DMFT calculation. The order parame- 
ters arc linear combinations of the number of particles, 
and can be easily computed from the Matsubara Green's 
functions obtained by solving the DMFT equations. For 
a given sublattice a and species a, 

^ oo 

nacr = ^ GQo-(«i^n) cxp(ia;„0"'") (Al) 

Taking into account the problem of convergence involv- 
ing \u)n\ — > OO, Eq. Al can be expressed in the more 
numerically tractable manner as 

11°° 

na<T = 2 + ^5 [Gaa{i(^n)] (A2) 



The potential energy can be simply computed as U 
times the double occupancy. In a general case away 
from half-filling, it might be also necessary to add a term 
= fJ^irih + ni). The specific heat is then numerically 
evaluated as 



^" dT 



(A6) 



The entropy requires the T integration of the specific 
heat, since at constant volume TdS — CydT. This can 
be achieved by integrating C^/T from T = or T — > oo. 
In general, the latter option is better since the asymp- 
totic value S{T — > oo) = A^ln(4) is already known, and 
therefore the expression becomes 



S{T) ^ 5(oo) 



dT 



,Cy{T') 



T' 



(A7) 



The thermodynamical observables, however, require 
the calculation of the internal energy E = {H). In the 
disordered phase, the kinetic term is given by 



(A3) 



Within DMFT for the Bethe lattice it can be rewritten 
as , 



2 giw„0+ 



(A4) 



However, in some cases, such as an ordered state, 
where there is a singularity in the specific heat, it is also 
necessary to carry out the integration from T = 0, which 
results in 



S{T) = S{Q)+ r dT'^''^^'^ 



T' 



(A8) 



To use this expression one needs a good justification for 
the value of S'(O). 
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